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We present measurements of various geometrical characteristics of monopole clusters in 5(7(2) 
lattice gauge theory. The maximal Abelian projection is employed and both infinite, or percolating 
cluster and finite clusters are considered. In particular, we observe scaling for average length of 
segments of the percolating cluster between self-crossings, correlators of vacuum monopole currents, 
angular correlation between links along trajectories. Short clusters are random walks and their 
spectrum in length corresponds to free particles. At the hadronic scale, on the other hand, the 
monopole trajectories are no longer random walks. Moreover, we argue that the data on the density 
of finite clusters suggest that there are long-range correlations between finite clusters which can be 
understood as association of the clusters with two-dimensional surfaces, whose area scales. 



I. INTRODUCTION 

The interest in monopoles in non- Abelian gauge theo- 
ries is mostly due to the dual superconductor mechanism 
of the confinement, for review see Q. The mechanism 
assumes condensation of magnetic monopoles in the vac- 
uum. The idea is supported by the lattice data and the 
phenomenology of the lattice monopoles is quite rich. Es- 
pecially, in case of SU (2) gluodynamics which we will 
concentrate on in this paper. 

Detailed theoretical interpretations of the data appear, 
however, difficult because monopoles are defined not di- 
rectly in terms of the original SU(2) fields but rather 
in terms of projected fields. The use of a projection is 
rooted in the fact that monopoles are intrinsically U(l) 
objects and there are infinitely many ways to select a 
U(l) subgroup for the monopole definition. In particu- 
lar, monopoles of the maximal Abelian projection (MAP) 
are defined in the following sequence of steps (which we 
describe in a somewhat simplified way). First, one finds 
the maximal Abelian gauge maximizing the functional: 



RAbel = 



(AI( X )) 2 + (AI( X )Y 



(1) 



where A"(x) are components of the gauge field. Finally, 
monopoles are defined as singularities of the A 3 {x) fields 
which violate the Bianchi identities: 



properties of the monopole clusters. In particular, it is 
important to distinguish between the percolating clus- 
ter and finite clusters. The percolating cluster fills in 
the whole of the lattice and is in a single copy for each 
field configuration. Finite clusters are characterized, in 
particular, by their spectrum as function of the length I. 
Commonly, one introduces the corresponding densities of 
the percolating and finite clusters: 



4 Pperc^ ^sites: I fin — 4 Pfin^ ^ site 



(3) 



where l pe rc,lfin is the total length of the corresponding 
clusters, N S i tes is the number of lattice sites and a is the 
lattice spacing. 

An important question is whether position of the 
monopoles is of any physical significance. On one hand, 
monopole trajectories are distinguished by singularities 
of the projected fields, see (01, and they are point like 
in terms of the projected fields. On the other hand, one 
can suspect that these singularities are artifacts of the 
projection. 

Probably, a priori the latter possibility looks more rea- 
sonable since the monopole definition involves a non-local 
gauge fixation, see However, there are accumulating 
lattice data which indicate that in the maximal Abelian 
gauge the monopoles, which by construction have size 
of the lattice spacing a, are physical. In particular, the 
density of the percolating monopoles scales: 



j™ n (x) = e^»°d v d p Al{x) 



(2) 



where j™ on (x) is the monopole current and all the ex- 
pressions are well defined on the lattice. 

Geometrically, monopoles are represented by closed 
trajectories (on the dual lattice) and one usually discusses 
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Pperc = 0.65(2) a 



3/2 
SU(2) 



(4) 



where crsu(2) is the string tension [21j and we quoted the 
data from Ref. Jsj where references to earlier papers can 
also be found. 

Let us emphasize that the result Q implies that the 
corresponding monopole trajectories are meaningful even 
on the scale a. Indeed, the total length of the percolating 
clusters which scales in the physical units is added up 
from small steps of size a and finally is independent on 
a. Thus, one could argue, see, e.g., Q, that there are 
gauge invariant objects which are detected through the 



2 



projection. It is of course an intriguing hypothesis which 
is worth to be thoroughly checked. Unfortunately, this 
suggestion cannot be scrutinized theoretically since the 
anatomy of the lattice monopoles in terms of the non- 
Abelian fields is largely unknown, for review see, e.g., 
Rather one should apply further phenomenological 
tests. 

In particular, one can measure non-Abelian action 
associated with the monopoles, see |j| and references 
therein. For monopoles in the percolating clusters it 
turns out that, at least at presently available lattices, the 
action corresponds to a monopole mass which diverges in 
the ultraviolet: 



M{a) 



const 



(5) 



M(a) is the lattice analog of the 'magnetic mass' in the 
continuum, M(a) ~ J H 2 <i 3 r where H is the mag- 
netic field. The result JSJ strongly suggests that the lat- 
tice monopoles are associated with singular non-Abelian 
fields and the singularity J5J is not an artifact of the pro- 
jection. 

The next question is how singular fields can be impor- 
tant at all since their contribution in the limit a — > is 
suppressed by an infinite action |22| . An apparent answer 
to this question is that the action and entropy factor are 
both divergent in the ultraviolet but cancel each other 
in such a way that, say, the total density of monopoles 
in the percolating cluster scales in the physical units, see 
(JIJ. An example of such a fine tuning is provided by 
U(l) theory, see, in particular, 0. In that case the fine 
tuning is ensured by choosing tuned values of the electric 
charge. In case of non-Abelian theories the hypothesis 
on the fine tuning 0, 13 is rather phenomenological and 
the mechanism of the fine tuning is to be understood yet. 

To summarize, 

(1) there is evidence that at presently available lattices 
the monopoles defined within the maximal Abelian pro- 
jection are associated in fact with singular non-Abelian 
fields fine tuned to the corresponding entropy factors. 

(2) This means that the structure of the non- 
perturbative vacuum fluctuations is quite different from 
the standard picture of 'bulky' fields of the size Aq ~ d . 

(3) Although the conclusion on the fine tuning is 
strongly suggested by the data like (@J| , (J5J further checks 
of it are highly desirable. Indeed, the whole basis is the 
phenomenology of the lattice monopoles and further data 
can bring further insight. 

Motivated by these considerations, we have inves- 
tigated in more detail the scaling properties of the 
monopole clusters 23] . The outline of the paper is as 
follows. In Sect. |n] we present detailed studies of the ge- 
ometrical elements of the percolating cluster. Namely, 
the cluster consists of self-crossings and segments be- 
tween the crossings. We report on the measurements 
of distribution in the length of the segments, long-range 
correlations between directions of the links occupied by 
monopole currents, scaling properties of the segments 




FIG. 1: The distribution N(l)/N, (3 = 2.4, the lattice size is 
24 4 . 



FIG. 2: N(l)/N for trajectories which are not closed loops; /3 
value and the lattice size are the same as in FigQ 




characteristics. In Sect. IIIII we first present data on 
the spectrum of the finite size clusters fsubsection llll Aj) 
and then on their density (subsection IIII B ll . In subsec- 
tion IIII CI we present results on the correlations of the 
monopole trajectories at large distances which are sen- 
sitive, at least in principle, to the glueball masses. Dis- 
cussions of the observations made in this work are in 
Sect. IIVI Finally, some details of simulations are given 
in the Appendix. 



II. SEGMENTS OF THE PERCOLATING 
CLUSTER 

A. Distribution in the segments length 

The percolating cluster consists of segments (that is, 
trajectories between crossings) and crossings. The seg- 
ments, in turn, are made from the links on the dual lat- 
tice. We will consider the length of the segment, l se gm 
and the Euclidean distance between the end points of the 
segments, d. 

In Fig.^we show the normalized distribution, N(l)/N, 
of the lengths of the segments; N(l) is the number of the 
segments with the length I in lattice units, N is the total 
number of segments. 

The oscillations of N(l)/N are clearly seen, the num- 
ber of trajectories with the even I is systematically larger 
than that with the odd I. This effect is due to the large 
number of the closed trajectories (loops), which are seg- 
ments with d = and I ^ 0. In Fig. |2|we show N(l)/N 
only for not closed trajectories, d 0. 
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FIG. 3: Correlations of the directions of the links lying on the 
monopole current; f3 = 2.6, the lattice is 28 4 . 



TABLE I: The fitted masses from data for P S ame 
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with Fig. □ 




FIG. 4: Average segment length, {l segm ), vs. lattice spacing 



B. Long-range correlation between links directions 

Next, we turn to angular correlation between the links 
on the monopole trajectory. Let us call link Co belonging 
to the percolating cluster as the "initial" one. Then one 
can measure the probability that the link C;, connected 
to the link Co by the monopole trajectory of length I, 
has the same direction. Of course this probability is a 
decreasing function of I, but it occurs that there exists 
a "long memory" of the initial direction. This fact is 
illustrated in Fig. |3| where the correlation of the direction 
of the initial link Co with the direction of the link Ci is 
shown. The direction of C; can be the same as that of Co, 
or opposite, or else (neither the same nor the opposite). 
We normalize all three probabilities in such a way that 
they are equal to unity for random walk if I ^= (for 
I = the opposite direction is forbidden). From Fig. it 
is seen that even for links separated by 17 lattice steps 
the most probable direction is the same as the direction 
of the initial link. For smaller values of j3 the correlations 
of the directions exist for not so large values of I. The 
deviation from unity of the probability to have the same 
direction falls off exponentially: P sam e = 1 + A s e~ fJ ' sla . 
The probability to have the opposite direction behaves 
as follows: P opP osit = 1 — A e~^° la . These fits are shown 
by dashed lines in Fig. El The masses [i s and fi coincide 
within the error bars and are independent of the lattice 
spacing, see Table [I] 



C. Scaling properties of the segments 

It was already mentioned in the Introduction that the 
density of the percolating cluster scales, that is p verc 
has a well defined value in the continuum. Here we will 
present further data on the scaling properties of geomet- 
rical elements of the percolating cluster. 

First, we observe that the average length of the seg- 
ment of the monopole trajectory between crossings, 
(lsegm), does not depend on the lattice spacing, see Fig.U 
Making fit by a constant for (3 > 2.35 we obtain 

(l segm ) = 1.60(1) fm. 

The way used to convert lattice results into physical units 
is explained in Appendix. The lattice spacing values a 
for various /3's are also given in Appendix, see Table III. 

In the average Euclidean distance between the cross- 
ings (d) the violations of the scaling are more significant, 
see Fig. 5. However, the deviations from the scaling 
can be approximated for /3 > 2.35 by linear in a cor- 
rections. Then in the continuum limit a — > 0, (d) has a 
non- vanishing value (~ 0.20 fm), as it is seen from Fig.0 

In Fig.|S|the number of crossing points per unit length 
of the monopole trajectory is shown by open circles. This 
number weakly depends on a and seems to have the con- 
tinuum limit (a — > 0) « 0.3 /m _1 . 

The situation changes if we exclude the closed loops of 
finite lattice length connected to the percolating cluster. 
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FIG. 5: (d) vs. lattice spacing. The dashed line shows the 
linear fit for j3 > 2.35. 



FIG. 7: The length distribution of finite clusters. (3 = 2.4, 
the size of the lattice is 32 4 . 
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FIG. 6: Number of crossings per unit physical length vs. lat- 
tice spacing in full percolating cluster and in the percolating 
cluster with excluded closed loops. Numbers 4, 6, ...12 in the 
legend mean that closed loops of the length up to 4,6, ...12 
are excluded. The dashed lines show the linear fits. 



The number of crossings reduces and is compatible with 
zero in the continuum limit if we exclude the loops of 
the length up to 8 (or up to larger length). This fact is 
illustrated in Fig. 

The continuum limit value (N cross ) /l perc = 0.3 fm 
for the full percolating cluster together with the data for 
the percolating cluster density (Fig. [HJ corresponds to ap- 
proximately 10 crossings in percolating cluster per hyper- 
cube 1 fm 4 . To complete the picture of the structure of 
the percolating monopole cluster in the continuum limit 
the data for (l S egm), (d) and (N cross /l perc ) can be com- 
pared with the average monopole radius (p m ) ~ 0.05 fm 
and average inter-monopole distance (R) ~ 0.5 fm [5j. 



III. FINITE CLUSTERS 
A. Length spectrum of finite clusters 

In Fig. [7| we show the number of the finite clusters of 
the length I (in the lattice units) vs. I. The dashed line in 
this figure is the fit of the data by function Const/l a . It 
is important that l/l 3 behavior seems to be valid for all 
considered values of the lattice spacings as it can be seen 
from TableHTIwhere we show results of the fit for various (3 
values. As we discuss in Section llV Cl the dependence 1 /I 3 
means that finite clusters develop in the four-dimensional 
space. 



B. Monopole densities 

As we discussed in the introduction the magnetic 
monopole clusters fall into two different classes yl ^| : 
" small" , or finite clusters which have finite size in lat- 
tice units, and "large" clusters which percolate through 
the lattice. It was demonstrated @ that the percolat- 
ing cluster is responsible for the string tension. If the 
size of the lattice is large enough, there is a natural dis- 
tinction between "small" and "large" clusters, since in 
each configuration there exists only one cluster which is 
much longer than all others and there is a gap in the clus- 
ter length distribution 0. If the lattice size is not large 
enough the percolating cluster decays into one or more 
large clusters plus several monopole currents which wind 
through the boundaries of the lattice . It occurs that 
the monopole density for the sum of these large clusters 
and the winding trajectories scales @- We simply call 
this collection of the monopole currents as percolating 
cluster. 

In Fig. [S] we show the monopole density p perc as a 
function of the lattice spacing a. Our results for p perc 
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TABLE II: Length spectrum of finite clusters, power fit parameter 





FIG. 8: Density of the finite clusters pfi„ and percolating 
clusters p per c', solid and dashed lines, are fits by a constant 
and Eq.0 respectively. 



obtained with higher statistics than in [2] agree within 
error bars with values obtained in [2| for coinciding values 
of [3 on a is rather weak and we fit it by a constant 
for j3 > 2.35 The resulting density of the percolating 
monopoles in the continuum limit is 



Pperc 



7.70(8) fm 



-3 



(6) 



This value is shown by the solid line in Fig. 03 It is in 
agreement with value obtained in |2| within error bars. 

The density of the finite clusters, on the other hand, 
is divergent for a — > and it can be fitted (dashed curve 
on FiglSJ as 



_ n C 2 

Pfin — Oi + — — 



(7) 



where 



Ci = -6.1(5) fm- A , C 2 = 1.55(4) fmr 2 . 

The negative value of the constant C± means that the fit 
is not valid for large values of the lattice spacing. 

We made also the local fits of pfi n by function A ■ a~ a , 
where A and a are fit parameters. The fit results for a are 
shown in Fig. To get a point on this figure used for the 
fit three neighbor data points in Fig. [5] Corresponding 
fit intervals are depicted in Fig.[5]as X axes error bars. It 
can be seen from Fig.[!|lthat the values of a are changing 
from the value close to 3 down to the value approximately 



FIG. 9: The dependence a(a). 



equal 1, i.e. at small a this fit becomes consistent with 
the fit Eq.Q. 

Another possible fit, suggested in Ref. |3j, pfi n = A- 
a~ a gives A and a strongly dependent on a, see Fig. [5] 
However, at smallest a we have the same 1/a behavior of 
the density. 



C. Correlator of the monopole currents 

We turn now to discussion of correlation functions of 
the monopole currents. We will consider three defini- 
tions: 

• Gi (x — y) is the probability that point x and point 
y are connected by a monopole trajectory. 

• G*2 (x — y) is the probability that point x and point 
y are connected by a monopole line belonging to 
the percolating cluster. 

• G3 (x — y) is the probability that the monopole cur- 
rent crosses point x and point y. 

Since the small clusters have finite size in the lattice 
units (see Sect. lIII Al the functions Gx(x — y) and G2{x — 
y) coincide with each other for large d = \x — y\. The 
asymptotic behavior for d = \x — y\ — > 00 is similar for all 
correlation functions and was suggested first in H3 (see 
also [13): 



G k (x-y) -> C k + A k exp{-m k d} 



(8) 
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FIG. 10: Polar plots for G 3 (d) at /3 = 2.3, the lattice size is 
16 4 . 



Finally, we determined the masses irik in Eq. ©. The 
fits made for the distance d > 0.3 fm show that the 
masses rrik do have scaling behavior, but we do not have 
the high precision in their determination: 

mi,2 = 1.87(16) GeF, m 3 = 2.4(5) GeV . 

Due to large errors we can claim only that the results do 
not contradict the prediction that rrik coincide with the 
mass of the lightest SU(2) ++ glueball as measured, in 
particular, in Refs. (13j |. 
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FIG. 11: The scaling behavior of the correlation function 
Gi(d). The constant Ci has been subtracted. 



IV. PERCOLATION LAWS AT SHORT AND 
LARGE DISTANCES 

In preceding sections we presented data on the phe- 
nomenology of the monopole clusters and have chosen 
to postpone comments on their implications till this sec- 
tion. It might worth emphasizing from the very begin- 
ning that not every piece of information that we got has 
a straightforward interpretation. Moreover, it is useful 
to distinguish between the properties of the clusters at 
relatively short distances, ^qqd ' a wnere ' is the 
length along the monopole trajectory, and at large dis- 
tances / > ^qcd- As we shall see, it is the short distance 
behavior which is quite well understood while for larger 
distances one could rather speak of accumulation of the 
data. This might look paradoxical since traditionally the 
non-perturbative fluctuations are deprived of interesting 
short-distance behavior. 



where k = 1,2,3 and rrik are expected to be the mass of 
the lightest glueball. Moreover, since by definition (see 
(PI) the monopole density is the probability to find the 
monopole current on the given link, the constants Ck are 
equal to the squares of the corresponding densities. It is 
obvious that Ci,2 = (pperc) and C3 = (p per c + Pfin) ■ 

The results of the measurements are presented in 
Figs. [HI In Fig- HOI we show the results of the ro- 
tational invariance check of the correlators Gk at large 
enough d. Namely, the dependence of G3 (x — y) on the 
angle between 4-vector x — y and one (arbitrary) axis in 
polar coordinates is shown. For illustrative purposes, we 
use different normalizations of G3 on different plots in 
Fig. ^3 The data show nice rotational invariance. Note 
that the gap in the angles seen for large d is due to the 
geometry of the lattice. 

One of the most important properties of the correla- 
tors is their scaling. The data do indicate that the cor- 
relators do have the scaling behavior. This is illustrated 
in Fig. 1111 where we show Gi(d) in physical units for 
various f3. Note that the dimension of p is fm~ 3 . Re- 
spectively, the dimension of the correlation functions Gk 
is fm~ 6 and they are sensitive to possible scaling viola- 
tions. Comparing the data for various volumes at (3 = 2.4 
we found the lattice volume independence of the correla- 
tors Gfc(d). 



A. Monopoles as free particles at short distances 

One of our main results is observation of scaling for a 
few observables : 

• average distance between the self-crossings of the 
percolating cluster, as measured along the trajec- 
tory, scales (see Fig. 0}; 

• average Euclidean distance between the same 
points approximately scales, see Fig. [5J 

• correlation length of direction alignment of links 
along the trajectory of the percolating cluster 
scales, see Table H] 

• correlators of the vacuum monopole currents scale, 
as illustrated in Fig. 1111 

• we confirm scaling of the density of the percolating 
cluster, see Fig. |S] 

Each of these checks is highly non-trivial since the most 
natural explanation of any scaling behavior is that the 
monopoles of size a are physically meaningful. Thus, 
each example of the scaling behavior confirms, at least 
indirectly, fine tuning of the lattice monopoles. 
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Now, let us discuss briefly the geometry of the 
monopole cluster at various distances. At the very 
short distances, of order a, monopole trajectories can be 
thought of as (approximate) random walk. At these dis- 
tances both the action and the entropy are of order 1/a. 
Let us illustrate this statement by a well known equa- 
tion (see, e.g., |l4() for the propagating mass of a scalar 
particle m; 



prop " 



(M(«) - f) 



(9) 



M(a) is the field theoretical mass of the same particle, 
discussed in the Introduction on the monopole example. 
A more general definition of M(a) is that the action of 
the particle is 

S d = M(a)-l , 

where I is the length of a path connecting two points be- 
tween which the particle propagates. Moreover, the In 7 
factor is due to the entropy. Eq. @ is derived in the ap- 
proximation of neglecting neighboring trajectories. Note 
a dramatic difference between M{a) and m pr0 p. Namely, 
to ensure a finite propagating mass m prop one needs an 
ultraviolet divergent mass M(a) (that is, the mass which 
determines the classical action). 

Upon accounting for possible cancellations between the 
action and entropy one can introduce an 'effective action', 



5', 



I ■ u , (/ > a) , 



which governs the behavior of the trajectories at dis- 
tances much larger than a. The very fact that the density 
of the percolating cluster scales implies that at distances 
of order a the divergent factors in the action and entropy 
cancel each other and fj, ~ ^-qcd- This is just the fine 
tuning hypothesis which we heavily rely on. 

Eq. © is useful for orientation in mass scales which 
came out from various fits to the data. In particular fits 
to the correlation of the links resulted in the attenuation 
factor of order exp(— /j, s l) where fi s is about 300 MeV 
and I is measured along the trajectory, see Table [I] At 
first sight, 300 MeV might look too low a mass. How- 
ever, if we apply for estimates Eq. JjJJ we would obtain 
rather m prop ~ (300 MeV) /a. In other words, we 
do have evidence that the most singular part of the ac- 
tion (as measured on the lattice, see the Introduction) 
is indeed cancelled by the entropy factor. But the value 
^ s w 300 M eV does not necessarily corresponds directly 
to a physical mass. 

Once this cancellation is confirmed by the data, one 
can make predictions about short clusters with the length 
^■qcd Namely, at such length one can 

neglect the mass factor and the spectrum of the closed 
loops should correspond to free particles. Then one can 
derive (see [j| and references therein): 



N(l) 



const 



R(l) ~V~l , 



where R(l) is the radius of the cluster of the length I. Let 
us emphasize that l|10|l is a consequence of free field the- 
ory in Euclidean space-time (actually, in d = 4 Eq. I|10|) 
survives adding Coulomb like interaction as well 8J ) . Eq. 
(|10fl is derived most naturally in the so called polymer 
representation, see, e.g., [sLU^I]. 

Measurements of N(l) and R(l) were reported in Q 
and are in perfect agreement with (|1UI) . Our analysis 
confirms (llOfl on larger statistics and for smaller values 
of a. In particular, the data on the I distribution are 
summarized in the Table[n] We did not include the data 
on the R(l) which we have but they do confirm (|10fl for 
all a tested. 



B. Percolating cluster 

Existence of a single infinite cluster is typical for all 
percolating systems in the supercritical phase, see, e.g., 
[lfij . There is a general theorem on the uniqueness of 
the percolating cluster at p > p c where p is the prob- 
ability to have a link 'open' (in our case, to belong to 
a monopole trajectory) and p c is the point of the phase 
transition. The percolating cluster is characterized by 
the probability 9(p) of a given link to belong to the clus- 
ter (in our case this probability is equal to the cluster 
density, p per c)- Generically, 9(p) ~ (p — p^ 13 , > 0. In 
our case, 



6(p) ~ (a ■ Aqcd) 3 « 1 



(11) 



Smallness 9(p) implies closeness to the point of phase 
transition to the percolation. In the limit a — > we hit 
exactly the point of phase transition. Expression Hll|) 
can be considered as formulation of the hypothesis on 
fine tuning of the monopoles in terms of the percolation 
theory 

Another general property of the percolating cluster is 
that its fractal dimension coincides with the dimension 
of the space: 



d p ;; c = d 



(12) 



Which is trivially satisfied in our case, since the length 
of the percolating cluster is proportional to the whole 
volume, V4. 

However, to the best of our knowledge, the percolation 
theory is not powerful enough to predict or explain the 
data which we have on such characteristics as l segm , d. 

It is worth emphasizing that geometrical characteris- 
tics of the percolating cluster at length of order {l aegm } 
differ quite drastically from the corresponding character- 
istics of the finite clusters. First of all, the monopole 
trajectories are no longer random walks. Indeed, for a 
random walk we would have had: 



(10) 



1 



(13) 
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where I is the distance between two points measured 
along the trajectory and d is Euclidean distance between 
the same points. 

Our results on l segm and d are given in Fig. 01 and 
respectively. The data certainly rule out l|13|) for any 
finite m pr0 p- Moreover if we assume that it is only the 
leading power of 1/a which is cancelled in Eq. @ then 

mriprop ~ y/A-QCD/a- Our data are not consistent with 
(|13fl for such mass either. 

Also, existence of correlation between directions of the 
links, (see Sect. IIIBf) . could not be reconciled with the 
random walk. 



C. Long-range correlation of finite clusters 

Turn now to discussion of the data on the density of 
the short clusters, see Sect. IIII Bl Expectations for pji n 
are easy to formulate. Indeed, the spectrum N(l) ~ Z~ 3 
implies that the average length of the finite clusters is 
saturated in the ultraviolet, i.e. at length of order a. 
Then the density of the finite clusters should not depend 
on Aqcd at all and on the dimensional grounds pfi„ ~ 
a" 3 . 

Instead, we find that the most singular piece in the den- 
sity pfi n ~ Aqcd ' a ■ The only interpretation of this 
striking observation is that even short, or finite clusters 
are not independent at large scales. Again, on dimen- 
sional grounds alone it is clear that the finite clusters are 
in fact associated with two-dimensional surfaces whose 
area scales. 

At first sight, this conclusion looks too bizarre. How- 
ever, in fact it is known to be true from independent 
measurements on the P-vortices. Namely, it was observed 
first in Ref. that monopoles are associated predom- 
inantly with P-vortices, whose total area scales. This 
conclusion was recently reinforced by measurements for 
smaller lattice spacings |l8| . 

Thus, we are invited to think about percolation of 
monopoles as a two-step process: monopoles percolate 
on the surface of dimension d = 2 and the surface itself 
percolates in d — 4 space-time. As far as we know, this 
type of percolation has not been studied theoretically at 
all. And we can, therefore, suggest only very preliminary 
considerations. 

First the probability 6(j>), see Eq. ff 1 II) is to be thought 
of as a product of two factors: 

d(p) ~ (a ■ Aqcd) 3 ~ (a-A QCD )-^ , (14) 

where the first factor is to be interpreted as a power of 
(p( 2 ) — pc 2 ') for percolation of monopoles on a surface 
while the second factor accounts for suppression of the 
'phase space' available on a d — 2 space, V2, as compared 
to the whole lattice volume V4. Note that for probabili- 
ties defined on the surface we still have according to i|14|) 

{ P {2) -pf ] ) ~ (a -Aqcd) « 1, a^O. 



Another major change, is that finite clusters now oc- 
cupy a finite fraction of the links on the surface. Further- 
more, one could also speculate that for finite monopole 
clusters of length I ~ Aqq D we would still have 

R ~ jV«i _ VI j d = 2 

As the last remark, let us notice that the interplay 
between the monopole trajectories and P-vortices is more 
complicated than simply saying 'monopoles belong to P- 
vortices'. Indeed, for short clusters the length spectrum 
is sensitive to the dimension of the space [3, Ha] : 



*(0 ~ jim > ( 15 ) 

and the spectrum ~ l~ 3 corresponds to d = 4. In this 
sense, at short distances the monopole trajectories are 
'primary' objects. Thus, one can expect that for larger 
I where the infrared behavior of the clusters sets in the 
spectrum is changed into l/l 2 . 



V. CONCLUSIONS 

There is accumulated evidence that monopoles defined 
within the Maximal Abelian projection appear physical 
even if studied with resolution of the lattice spacing a. In 
particular, we have seen that there are amusingly simple 
scaling laws for various observables which depend only 
on the product (a ■ Aqcd) and look perfectly gauge- 
invariant. Thus, non-perturbative fluctuations in the lat- 
tice SU (2) appear to have non-trivial structure in the 
ultraviolet as well as in the infrared. 

What is specific for the field, is that phenomenology 
seems to be far ahead of theory. In particular, we have 
argued that there are indications that percolation of the 
monopoles is to be considered as combination of perco- 
lation of the trajectories over a surface d = 2 and of per- 
colation of the surface over the d = 4 space. To the best 
of our knowledge, no theoretical framework is available 
to describe such a percolation. 
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2.30 2.35 2.40 2.40 2.40 2.40 2.45 2.50 2.55 2.60 

a(fm) 0.1655(13) 0.1394(8) 0.1193(9) 0.1193(9) 0.1193(9) 0.1193(9) 0.0996(22) 0.0854(4) 0.0713(3) 0.0601(3) 

L 16 16 16 24 28 32 24 24 28 28 

# of conf. 100 100 300 137 14 35 20 50 40 50 



TABLE III: SU(2) configurations 



Appendix 

The list of configurations used in this work is given 
in Table IrTTl To fix the Maximal Abelian gauge we use 
the Simulating Annealing algorithm 0] and study 10 
randomly generated gauge copies for each configuration 
to avoid the Gribov copy problem. 



Most of the quantities in the paper are given in physical 
units (e.g. in fm). In order to express the lattice data 
in physical units we use the data of Refs. [2£j for the 
SU(2) string tension and assume that y/a = AAQMeV. 
The corresponding values of the lattice spacing are given 
in Table ITTT1 



[1] M.N. Chernodub, M.I. Polikarpov, in 'Cambridge 1997, 
Confinement, duality, and nonperturbative aspects of 
QCD \ p. 387; hep-th/9710205 

[2] V. Bornyakov and M. Miiller-Preussker, Nucl. Phys. 
B106 (2002) 646. 

[3] A. Hart and M. Teper, Phys. Rev. D58 (1998) 014504; 
Phys. Rev. D60 (1999) 114506. 

[4] S. Coleman, in "The Unity of the Fundamental Inter- 
actions", Erice lectures, 1981, ed. A. Zichichi, Plenum, 
London (1983) p. 21; 

M.N. Chernodub, F.V. Gubarev, M.I. Polikarpov, V.I. 
Zakharov, Phys. Atom. Nucl. 64 (2001) 561, Yad. Fiz. 
64 (2001) 615, (hep-th/0007135|. 
[5] B.L.G. Bakker, M.N. Chernodub and M.I. Polikarpov, 
Phys. Rev. Lett. 80 (1998) 30; 

V.G. Bornyakov, F.V. Gubarev, M.I. Polikarpov, 
T. Suzuki, A.I. Veselov and V.I. Zakharov, Phys. Lett. 
B537 (2002) 291. 
[6] A.M. Polyakov, Phys. Lett. B59 (1975) 82; 

T. Banks, R. Meerson, J. Kogut, Nucl. Phys. B129 
(1977) 493. 

[7] V.I. Zakharov, "Hidden mass hierarchy in QCD", 
hep-ph/020404. 

[8] M.N. Chernodub, V.I. Zakharov, "Towards un- 
derstanding structure of the monopoles clusters", 
hep-th/0211267 

[9] P.Yu. Boyko, M.I. Polikarpov and V.I. Zakharov, "Geom- 
etry of percolating monopole clusters", hep-lat/0209075 
[10] S. Kitahara, Y. Matsubara and T. Suzuki, Progr. Theor. 
Phys. 93, 1 (1995); 

M. Fukushima, A. Tanaka, S. Sasaki, H. Suganuma, 
H. Toki and D. Diakonov, Nucl. Phys. Proc. Suppl. 53, 
494 (1997). 

[11] T. L. Ivanenko, M. I. Polikarpov and A. V. Pochinsky, 
JETP Lett. 53 (1991) 543 [Pisma Zh. Eksp. Teor. Fiz. 
53 (1991) 517]; 

T. L. Ivanenko, A. V. Pochinsky and M. I. Polikarpov, 
Phys. Lett. B302 (1993) 458. 
[12] K. Langfeld, H. Reinhardt, "Monopole - anti-monopole 
excitation in MAG projected SU(2) lattice gauge theory", 
hep-lat/0206021 



[13] C. Michael and M. Teper, Nucl. Phys. B305 (1988) 453; 

C. Michael and S. Perantonis, J. Phys. G18 (1992) 1725; 
S. Booth et al (UKQCD), Nucl. Phys. B394 (1993) 509. 

[14] J. Ambjorn, "Quantization of geometry", 

hep-th/9411179 
[15] A.M. Polyakov, "Gauge Fields and Strings", Harwood 

academic publishers, (1987). 
[16] C. Itzykson, J.-M. Drouffe, "Statistical Field Theory", 

vol. 1, Cambridge University Press, (1989); 

D. Stauffer, A. Aharony, "Introduction to percolation the- 
ory" London et al.: Taylor and Francis 1994; 

G. Grimmett, "Percolation" Berlin et al.: Springer 1999, 
Grundlehren der mathematischen Wissenschaften, Vol. 
321. 

[17] J. A mbjorn, J. Giedt, J. Greensite, JHEP 9903 (1999) 

019, hep-lat/9903023 
[18] F.V. Gubarev, A V. Kovalenko, M.I. Polikarpov, S.N. 

Syritsyn, V.I. Zakharov, "Fine tuned vortices in lattice 

SU(2) gluodynamics, hep-lat/0212003 
[19] G. S. Bali, V. Bornyakov, M. Muller-Preussker and 

K. Schilling, Phys. Rev. D54 (1996) 2863. 
[20] J. Fingberg, U. M. Heller and F. Karsch, Nucl. Phys. B 

392 (1993) 493. hep-lat/920812 

G. S. Bali et.al. I nt. J. Mod. Phys. C4 (1993) 1179-1193 
hep-lat/9308003 

G. S. Bali, K. Schilling, A. Wachter Phys. Rev. D 55 
(1997) 5309-5324 hep-lat/9611025 

G. S. Bali, K. Schilling, C. Schlichter Phys. Rev. D 51 

(1995) 5165-5198 hep-lat/9409005 

B. Lucini, M. Teper JHEP 0106 (2001) 050 

|hep- lat/010 3027l 

[21] In estimates we will also use Aqcd although we study 
actually SU(2) gluodynamics, o- S u(2) ~ ^qcd- 

[22] Hereafter we tacitly assume that the a dependence ob- 
served at the presently available lattices will persist at 
smaller lattice spacings as well. The smallest a available 
now is about a « (3 GeV)~ l . 

[23] The preliminary numerical results were presented at the 
"Lattice 02" conference |<J. 



